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Abstract 

We use elementary variational arguments to prove, and improve on, gap 
estimates which arise in simulating quantum circuits by adiabatic evolution. 

There are several models for quantum computation The quantum Turing 
machine model and the quantum circuit model, are equivalent in the sense that any 
algorithm that runs in polynomial time in one requires only polynomial time in the 
other. There are also several "one-way" measurement-based models jSj, such as the 
cluster state model ^j, which can simulate any polynomial time quantum circuit in 
polynomial time. 
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In 12], Farhi, et al introduced quantum computation by adiabatic evolution of a 
Hamiltonian and showed that it could be simulated by a suitable quantum circuit. 
In this model, the time required is assumed to depend inversely on the square of the 
eigenvalue gap.^ In PO, a method was given for simulating an arbitrary quantum 
circuit with L gates by the adiabatic evolution of a Hamiltonian in time which is 
polynomial in L, using a simple "clock" model to construct the Hamiltonian. Some 
modifications jlj [7] of the Hamiltonian construction have been considered without 
introducing different techniques for estimating the gap. 

In P, the techniques used to prove the gap estimates are rather complicated. In 
this note we show that very elementary techniques suffice, and that one of the bounds 
can be improved by a factor of L. To make precise statements, we need some notation. 

Let |efc) denote the standard basis for C^; in particular, |ei) = (1,0, . . .0)^. Let 
—Ad denote the discrete Laplacian with Neumann boundary conditions, i.e.. 
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The Hamiltonian used in Lemma 3.5 of [Tj can be written as 

Hois) = si-Ad) + (1 - s)Id - (1 - s)|ei)(ei| 
and the block diagonal Hamiltonian in Lemma 3.12 as 



j=0 



(2) 



(3) 



where 



H^is) = si-Ad) + (1 - s)Id + [b, - (1 - s)]|ei)(ei|. (4) 

^This assumption does not take into account the fact that higher order terms in the asymptotic 
expansion may be needed, or their possible growth as the Hamiltonian changes with the size of 
the problem |S]. Although this may affect the time estimates, it seems unlikely to do more than 
change the order of the polynomial. In this note, we deal only with gap estimates and not with time 
estimates. For further details and references on issues involved see the report of the workshop at 
http:// www.perimeterinstitute.ca/ activities / scientific /PI-WORK-6 / related Jinks. php] 
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with bj an integer > 1 for j > 1, and 60 = 0. Let Ai(s) < A2(s) . . . denote the 
eigenvalues of Ho{s) and Ak{s) the eigenvalues of EI(s) also in increasing order. One is 
interested in g{s) = A2(s) — Ai(s) and G{s) = ^2(3) — Ai(s), the energy gaps between 
the two lowest states of these Hamiltonians. 

The eigenvalue equation Hj\u) = X\u) written in terms of the vector components 
Uk of \u) is 

+ = (1 - A)Mfc k = 2,3,...d-l (5) 
subject to the boundary conditions 

+ f - A) = f M2 (1 - f - A)Md = (6) 

This is a second order difference equation with constant coefficients subject to bound- 
ary conditions. It can be solved exactly by elementary methods entirely analogous to 
those used to solve the "particle in a box" boundary value problem. For general s, bj 
the algebra can become somewhat tedious and we need only bounds on the lowest 
eigenvalues, not the full spectrum. Good estimates on the gaps can be obtained from 
one special case and a simple variational argument. 

Lemma 1. The lowest eigenvalue fiQ of —Ad + ||ei)(ei| satisfies 1 > yUo > ^• 

Theorem 2. The energy gaps for the Hamiltonians given by (j21) and Q are both 
0{d-^); mfact, g{s) > l/2d^, and G{s) > l/2d^. 



When d = L + 1, part (a) is equivalent to Lemma 3.5 of [T] and part (b) improves 
Lemma 3.12 of PP by a factor of L. Since we are interested in large d, we will henceforth 
not distinguish between estimates involving d and d ±1. 

Proof of Theorem [2J Consider two simple trial functions using the ground states 
of Hq at the endpoints s = 0, 1. First, 

(ei, Ho{s) ei) = s(ei, (-A^) ei) = |s. (7) 

Let \v) denote the normalized constant vector with elements "^fc = Then 

{v, H,{s)v) = {l-s)^ = {l-s)-^,{l-s). (8) 

Thus, Ai(s) < min{|s, (1 — s)^^} and the two curves cross at Sc where | < Sc = 
|^5f < I for c? > 2. Next, note that by the max-min principle 

A2(s) > ini {u,Ho{s)u) = {I- s) + sfXQ. (9) 
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Since /iq < |, one finds that |(2 — 3s) + s/zq is decreasing on [0, sj and attains its 
minimum at Sr- Thus 



A2(s) - Ai(s) > 



|(2 - 3s) + s/io s < Sc 

i(l - S) + s/io S > 



> 



> 



1 



(10) 



This proves part (a). To prove part (b) observe that the condition bj > 1 imphes that 
for j > 1 the Hamiltonian satisfies 

Hj{s) > s(-Arf) + (l-s)/d + s|ei)(ei| 

> s(-Arf + i|ei)(ei|) + (l-s)/d 

> s/io + (l-s). (11) 

Now, A2(s) is the minimum of A2(s) and the lowest eigenvalue of Hj{s) with j > 1. 
Therefore, © and (fTTj) imply that A2(s) > sfiQ + (1 — s) as well. Since Ai(s) = Ai(s), 
the argument above implies that gap for HI(s) satisfies G{s) > l/2d?. 



Proof of Lemma [2 The eigenvalue problem for — + ^|ei)(ei| is equivalent to 
setting s = 1, and 6 = | in (0) and We first look for solutions of the form 

Uk = e'^^ - e-'^^ with < A < 2. When A = 1 
condition are satisified. One can verify that for 9 — ^^^^ 
holds for n = 1, 2, . . . , ci. Since this gives d linearly independent solutions in the range 
< A < 2, the d eigenvalues are 1 — cos (|g)7r for n = 1, 2, . . . , (i. The smallest 
eigenvalue is 



cos^^, (0) and the first boundary 
^"~^7r, the second condition in (EI) 



/^o = l-cos(^). = ^-0(l)>l. 



for d sufficiently large. QED 



As a final remark, we note that — A^ = \X^X with X 



( 1 




-1 
1 -1 



'121 



\ 





so that X\vL) has elements — Uk+i and 



{u,-Adu) = UXu,Xu) 



d-l 



Uk+l\ 



k=l 



This lends itself to interpreting — A^ as a lattice analogue of the kinetic energy, rather 
than as a 3-local potential as in pp. Moreover, the interpolating Hamiltonians Hj 
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are linear combinations of —A, the Identity /, which can shift the spectrum but has 
no effect on the gap, and |ei)(ei| whose only effect is to modify the first boundary 
condition. 
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